DEM_TWI_Mississauga

library(fastmap)
library(raster)
library(ggplot2)
library(rgl)
library(plot3D)
library(rasterVis)
library(animation)

Elevation of Mississauga plotted from raster data

elevation <- raster("DEM.tif")
par(mar = c(4, 6, 4, 2) + 0.1)

plot(elevation, main = "Elevation of Mississauga (Meters)", xlab = "Longitude", ylab = "Latitude")

# DEM of Mississauga Plot Perspectives

persp(elevation, 
      phi=35, 
      xlab="Longitude", 
      ylab="Latitude", 
      zlab="Elevation", 
      ticktype = "detailed", 
      nticks = 3, 
      col = topo.colors(4), 
      border = NA, 
      shade= 0.4, 
      expand = 0.30)

persp(elevation, 
      phi=1, 
      xlab="Longitude", 
      ylab="Latitude", 
      zlab="Elevation", 
      ticktype = "detailed", 
      nticks = 3,
      col = topo.colors(4), 
      border = NA, 
      shade= 0.4, 
      expand = 0.30)

persp(elevation, 
      theta=45, 
      phi=1, 
      xlab="Longitude", 
      ylab="Latitude", 
      zlab="Elevation", 
      ticktype = "detailed", 
      nticks = 3, 
      col = topo.colors(4), 
      border = NA, 
      shade= 0.4, 
      expand = 0.30)

persp(elevation, 
      theta=90, 
      phi=1, 
      xlab="Longitude", 
      ylab="Latitude", 
      zlab="Elevation", 
      ticktype = "detailed", 
      nticks = 3, 
      col = topo.colors(4), 
      border = NA, 
      shade= 0.4, 
      expand = 0.30)

persp(elevation, 
      theta=135, 
      phi=1, 
      xlab="Longitude", 
      ylab="Latitude", 
      zlab="Elevation", 
      ticktype = "detailed", 
      nticks = 3, 
      col = topo.colors(4), 
      border = NA, 
      shade= 0.4, 
      expand = 0.30)

persp(elevation, 
      theta=180, 
      phi=1, 
      xlab="Longitude", 
      ylab="Latitude", 
      zlab="Elevation", 
      ticktype = "detailed", 
      nticks = 3, 
      col = topo.colors(4), 
      border = NA, 
      shade= 0.4, 
      expand = 0.30)

# Elevation Histogram (Unfilled sinks)

elevationunfilled <- raster("DEM.tif")

hist(elevationunfilled, main = "Distribution of Elevation vs Frequency (Unilled Sink DEM)", xlab = "Elevation", ylab = "Frequency")

Elevation Histogram (filled sinks)

elevationfilled <- raster("Fill_DEM.tif")

hist(elevationfilled, main = "Distribution of Elevation vs Frequency (Filled Sink DEM)", xlab = "Elevation", ylab = "Frequency")

Topographic Wetness Index Plot

TWI <- raster("TWI.tif")
par(mar = c(4, 6, 4, 2) + 0.1)


plot(TWI, main = "TWI of Mississauga", xlab = "Longitude", ylab = "Latitude"
)